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Quarkonia and B-Physics are among the first areas to be investigated with the first data collected by ATLAS. 
The ATLAS detector at CERN's LHC is preparing to take data from proton-proton collisions expected to start 
by the end of 2009. Investigation of the decay of B-hadrons represents a complementary approach to direct 
searches for physics beyond the Standard Model. Early B-Physics data will provide valuable information on 
the detector performance, as well as allow calibration studies in support of New Physics searches. Meaningful 
quarkonia studies performed with early data are expected to have the reach to make authoritative statements 
about the underlying production mechanism and provide cross-sections in this new energy regime. We review 
various aspects of prompt quarkonium production at the LHC: the accessible ranges in transverse momentum and 
pseudo-rapidity, spin alignment of vector states, separation of color octet and color singlet production mechanism 
and feasibility of observing radiative Xc decays. 



1. Introduction and motivations 

The main goal of ATLAS [T], one of the gen- 
eral purpose Large Hadron Collider (LHC) ex- 
periments, is looking for signals of New Physics 
beyond the Standard Model (BSM). While the ex- 
pected BSM signatures are mostly at high trans- 
verse momenta (pt), the B-Physics processes are 
concentrated in the low px region. Nevertheless 
ATLAS is able to provide B-Physics measure- 
ments based on its precise vertexing and tracking, 
its muon identification and a dedicated trigger 
strategy. The LHC is expected to produce a large 
number of B decays and heavy quarkonia states 
such as J and T in the low luminosity runs dur- 
ing the first few years of running. Their sizable 
branching fraction into charged lepton pairs al- 
lows for easy separation of these events from the 
expected huge hadronic background. The AT- 
LAS efforts concentrate on those B decays that 
can be selected by the first and second trigger 
levels. The most favorable B and heavy quarko- 
nium trigger signatures will be events decaying to 
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/x pair, either directly or via a J ftp meson. Us- 
ing the J/ij) trigger ATLAS will be able to col- 
lect unprecedentedly high statistics of B s — » Jtpcf) 
decays, allowing measurements of CP violating 
effects, which are predicted by some BSM mod- 
els to be significantly larger than the SM pre- 
diction. The di-muon trigger will also give AT- 
LAS access to potentially large numbers of rare 
B s —> fj,fx decays. Exclusive decay channels of 
B + and B® with two muons in the final state, 
will be measured for studying mass and lifetime 
as well as production cross-section. Being nar- 
row resonances, the quarkonia are perfectly suited 
for alignment and calibration of the ATLAS trig- 
ger and tracking systems. On top of that, un- 
derstanding the details of the prompt charmo- 
nium production is a challenging task and a good 
testbed for various QCD calculations, spanning 
both perturbative and non-perturbative regimes 
through the creation of heavy quarks in the hard 
process and their subsequent evolution into phys- 
ical bound states. 

2. ATLAS B-Physics di-muon trigger 

ATLAS uses a fast and efficient muon based 
B-Physics triggering strategy. The LHC will col- 
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lide two 7 TeV proton beams at a rate of 40 MHz, 
which together with pile-up will result in an inter- 
action rate of up to around 1 GHz. The ATLAS 
trigger acts to reduce this rate to around 100 Hz 
in three successive levels, whilst keeping only 
events 'of interest'. The trigger first level (level- 1) 
decision is based on coarse granularity of two sub- 
detector systems: the muon trigger chambers and 
the calorimeters. The second level trigger (level- 
2) receives from level- 1 data restricted to limited 
Regions of Interest (Rol). For a level- 1 muon, the 
level- 2 will use the information from the muon 
and inner detector tracking chambers to improve 
the muon momentum estimate, which allows a 
tighter selection based on this quantity. There are 
two specific types of trigger dedicated to quarko- 
nium: one which requires two level- 1 Rols cor- 
responding to two muon candidates with px val- 
ues above thresholds of 4 and 6 GeV respectively, 
and the other requires a single level- 1 Rol above 
a threshold of 4 GeV and searches for the second 
muon of opposite charge in a wide Rol at level-2. 
A complementary trigger on 3/ip> events is trig- 
gering on a single higher pt (10 GeV) muon and 
searching for a matched track at the off-line anal- 
ysis level. Before incorporating trigger and recon- 
struction efficiencies, the predicted cross-sections 
for pp —> Jjip —> + X were calculated for 

a number of px thresholds on the di-muon trig- 
ger (see Table Q] for details) . Figure Q] illustrates 
the distribution of cross-sections across the val- 
ues of the pt of the harder and softer muon from 
the quarkonium decay without any muon cuts ap- 
plied and zero polarization. The lines overlaid on 
the plots represent the muon trigger thresholds 
Pt > 6 GeV and pt > 4 GeV for the harder and 
the software muon respectively (denoted further 
by "6+4 GeV" or /i6/i4). Similarly "4+4 GeV" 
Gu4//4) refers to the trigger threshold pt > 4 GeV 
applied on both muons. The " 10+0.5 GeV" QitlO) 
refers to trigger pt threshold of 10 GeV applied 
on one of the muons only. In all cases the pseu- 
dorapidity, ?y, of a muon lies within an interval 
M < 2.5. 

Even if the bulk of the J/ip are produced un- 
der the trigger thresholds, nevertheless, due to 
the high cross-section the number of accepted J/ip 
will be substantially larger than at the Tevatron. 
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Figure 1. Density oiJ/tp production cross-section 
as a function of the hardest and softer muon pt 
of muons coming from 3 ftp. No cut was placed 
on the generated sample, but the overlaid lines 
represent the thresholds of observed events with 
trigger cuts applied. 



In the T case, due to it's relatively larger mass 
the bulk of the production is in the region near 
the muon threshold of 5 and 4 GeV. This means 
that by lowering the threshold for the higher mo- 
mentum muon we can significantly increase the 
number of recorded T events. 

Estimated quarkonia cross-sections for the 
three scenarios are presented in Table [T] The 
production will be dominated by J/ip and T(IS'). 
while tp' and T(2S), T(35) are expected to give 
smaller contributions. 

Precise cross-section measurements require a 
good understanding of the event and trigger selec- 
tion efficiency. To evaluate the trigger efficiency 
we use tag-and-probe (TAP) method applied on 
the real J/ip — > [i + fi" data. In the TAP method 
a single triggered muon from a reconstructed di- 
muon decay of an identified specific particle is 
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Table 1 

Predicted cross sections for various quarkonium 
states into muons, for different pt thresholds (be- 
fore trigger and reconstruction efficiencies). The 
last column shows the overlap between the di- 
muon and single muon samples. 
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(six times the expected average mass resolution) 
of the expected mass, the pair is considered as a 
charmonium candidate. 

The radial displacement of the two-track vertex 
from the beamline is used to distinguish between 
prompt J/ip originated at the proton proton inter- 
action and the B-hadron decays into quarkonia. 
The pseudo-proper decay time, r, is defined as 

_ L X y ■ Mj/^ 

T p T {J/i>) ' 1 } 

where Mj/^ and pr(J/ip) represent the J/ip in- 
variant mass and transverse momentum, and L xy 
is the transverse decay length of the meson. The 
resolution in the pseudo-proper decay time is ex- 
pected to vary from 0.110 ps for the low px char- 
monia down to 0.07 ps for the higher pt- As 
demonstrated in Figure [2] a cut on this quantity 
can efficiently distinguish between prompt and in- 
direct J /ip events. 

In the case of B + — > J/ipK + all tracks with 
positive charge which do not originate from the 
primary vertex are considered as K + . We fit the 
three tracks corresponding to a J/ip and the K + 
candidates to a common vertex. The momentum 
sum of the two has to point to the primary vertex. 
A maximum likelihood (ML) method is used to fit 
the B + mass distribution taking Gaussian prob- 
ability density function for the signal and linear 
distribution for the background. The mass resolu- 
tion in this channel is estimated, for an integrated 
luminosity of 10 pb- 1 to be cr = 42.2 ± 1.3 MeV 
with a total efficiency of 29.8 ± 0.8%. The differ- 
ential cross-section is obtained in this channel in 
different bins of px(B + ) using the number of the 
signal events and the efficiency from the ML fit to 
the invariant mass distribution. The total cross- 
section is determined with the same procedure 
but for the whole pr region of B + . For an inte- 
grated luminosity of 10 pb -1 , we determined that 
the total production cross section can be mea- 
sured with a statistical precision better than 5% 
and the differential cross-section can reach a pre- 
cision of 10%. The systematic uncertainties are 
expected to be dominated by the luminosity and 
the branching ratio. 

The channel B s — > J/ip<p has the potential to 
reveal New Physics signals. The weak mixing 



used as a tagged muon while we probe the trigger 
efficiency of the second muon. We verified that 
the efficiency derived with the TAP method on 
simulated events is in good agreement with the 
single muon trigger efficiency directly estimated 
from Monte Carlo samples. 

3. B-Physics early measurements 

The large production cross section of bb events 
will enable, even during the initial low luminos- 
ity phase of the LHC, the observation of vari- 
ous exclusive channels such as B + —> J/ipK + , 
Bd — > J/ipK° or B s — > J/ip(p. The measurement 
of the known properties of the B hadrons such as 
masses and lifetime will be used to study and tune 
the detector. The known channels will serve as a 
reference while searching for the rare B decays. 

The charmonium events and the B decays to 
J/ip share a similar detection strategy. In any 
event which passes the di-muon trigger, all the re- 
constructed muon candidates are combined into 
oppositely charged pairs, and each of these pairs 
is analyzed in turn. The invariant mass of the 
combined muons is required to be within a de- 
fined window around the J/ip mass. In the direct 
charmonium search, if the invariant mass of the 
two muons is above 1 GeV, we first attempt to 
refit the tracks to a common vertex. If a good 
vertex fit is achieved, the pair is accepted for fur- 
ther analysis. If the invariant mass of the refitted 
tracks is within 300 MeV for J/ip or 1 GeV for T 
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Figure 2. Upper - Pseudo-proper decay time, t, 
distribution for reconstructed prompt J/tp (cross- 
hatched, centred at zero) and the sum of prompt 
and indirect J /ip from B-decays (hatched, expo- 
nential distribution). 

Lower - Sources of low invariant mass di-muons, 
reconstructed with a /i6/i4 px trigger, with the 
requirement that both muons arc identified as 
coming from a primary vertex and with a pseudo- 
proper decay time smaller than 0.2 ps. The white 
area represents the background that was rejected 
with the pseudo-proper decay time treatment. 



phase 4> s , which in the SM is very small, may 
be enhanced by BSM processes. The time depen- 
dent angular distribution of this decay depends on 
seven physics parameters: two independent com- 
plex transversely amplitudes, the mean lifetime 
and the mass eigenstate width difference (Ar), 
and the weak mixing phase. The full analysis will 
require a fit to these parameters and is sensitive 
to statistics, experimental resolution of lifetime, 
mass and decay angles as well as flavour tagging 
performance and background rejection. The ex- 
perience of the Tevatron experiments taught us 
that a few thousands events only allow a two di- 
mensional profile likelihood fit in 4> S -AT plane. 
The LHC after the high luminosity runs will en- 
able a simultaneous determination of all seven 
parameters. With early data the ATLAS pro- 
gram will begin with calibration measurements 
supporting this analysis. The topological identi- 
cal channel Bd — > J/ipK°, which is expected to 
yield 15 times more statistics, is the the primary 
background but is also an essential as a control 
sample. It will allow high precision tests of the 
lifetime measurement as well as flavour tagging 
analysis calibration. With the earliest data, AT- 
LAS will use a fit to simultaneously access the 
mean mass and lifetime of the B s and the Bd 
mesons. This will enable the testing of our under- 
standing of the tracking at low px scale after 150 
and will start to improve the world preci- 
sions on these measurements after about 1 fb . 
Early on, topologically similar backgrounds will 
be admitted to the analysis and secondary vertex 
cuts will not be applied. After 10 the pre- 

cision on the Bd lifetime is expected to be 10% 
and similar precision for the B s mean lifetime is 
expected after 150 pb^ 1 . 

4. Charmonium analysis 

Charmonium production comprises of three 
main processes: direct singlet production, octet 
production and singlet production of x states. 
Each of these processes is characterized by dif- 
ferent differential cross-sections, 

The reconstruction efficiency of charmonium 
varies with px and r\. When the px(J/^p) > 
10 GeV we get a sharp rise in its acceptance. The 
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rise in the T case is less sharp and it reaches 
a high acceptance level around 30 GeV. Both 
channels reach a similar plateau at acceptance of 
around 85%. The trigger requires both muons to 
pass the pt threshold and so the angular separa- 
tion between the two muons is not large for most 
of the accepted events. We describe the open- 
ing angle by SR — y/6(f> 2 + Sif 2 , where 5<j> and Sn 
are the differences in the azimuthal angle and the 
pseudorapidity of the two muons from 3/tp. Typ- 
ical values for SR are ps 0.47 which means that re- 
quiring both muons to be above the pr threshold 
forces them to fly very close one to each other. 
For this reason the J/tp angular acceptance fol- 
lows closely the individual muon distribution and 
its dependence on material and detector effects. 
In the T case, due to its higher mass, it tends to 
be produced with higher px- Therefore the sep- 
aration angle SR is much wider, the muons do 
not go to the same area in the detector and the 
efficiency dependence on pseudorapidity is much 
smoother than in the J/tp case. The main sources 
of low invariant mass di-muons expected to domi- 
nate the background for prompt charmonium are: 
decays in flight of tt^ and K (which are reduced 
after cuts to below the 1% level); di-muon pro- 
duction via the Drell-Yan process (which are ren- 
dered essentially negligible by the trigger require- 
ments); continuum of muon pairs from beauty 
(and charm) decays and the main source - indi- 
rect J/ip production - the reduction of which is 
discussed below. 

All the background sources apart from the 
Drell-Yan pairs contain muons which do not orig- 
inate from the interaction point, and this is used 
to suppress their contamination by rejection of 
the events containing a secondary vertex if iden- 
tified. As detailed in the previous section a cut 
on the pseudo-proper decay time is used to dis- 
tinguish between the direct and indirect charmo- 
nium cases. Figure [5] shows the pseudo-proper 
decay time distribution of the two samples and 
the impact on the mass distribution when re- 
jecting events with a pseudo-proper time above 
0.2 ps. 

Using the /ilO trigger each reconstructed sin- 
gle muon candidate is combined with oppositely- 
charged tracks reconstructed in the same event 



within a cone of SR = 3.0. Any track, including 
those that were not identified as muons, are ex- 
amined. As in the di-muon analysis, we require 
that both the identified muon and the track are 
flagged as having come from the primary vertex. 
In addition, we impose a cut on the transverse im- 
pact parameter do, \do\ < 0.04 mm on the iden- 
tified muon and |c?o| < 0.10 mm on the second 
track, in order to further suppress the number 
of background pairs from B-decays. We obtain 
a invariant mass resolution close to that in 
those events passing the trigger requiring both 
muons to be confirmed. It's worth noting that the 
signal-to-background ratio around the J/"0 peak 
improves slightly with increasing transverse mo- 
mentum of J/ip- At higher pt the cos 9* accep- 
tance also becomes broader, which should help 
independent polarization measurements. 

ATLAS is capable of detailed checks of the pre- 
dictions of various models describing the evolu- 
tion of heavy quark antiquark pair into a quarko- 
nium bound state. Among the models which 
are considered are: Color Evaporation Model 
(CEM)0, the Color Singlet Model (CSM) [3] 
which failed to predict the J/ip production rate 
measured by CDF, and the Nonrelativistic QCD 
(NRQCD) Colour Octet Model (COM) g] which 
was proposed to explain this discrepancy. The 
good description of the Tevatron data by the 
COM model shown in Figure [3] is at least in part 
due to the tuning of some of its parameters which 
were determined from the same data. However, 
the difference between the CDF measurement of 
J/ip polarization dependence on transverse mo- 
mentum (pt) and the theory predictions (Fig- 
ure 0]) is motivating us to repeat the measure- 
ment in the LHC's higher energy and luminosity 
regime. 

The polarization of quarkonia can be measured 
using the angular distribution of the daughter 
particles produced in the decay. 



where cos 9* is the angle between the direction of 
the positive (by convention) muon from quarko- 
nium decay in the quarkonium rest frame, and 
the direction of the quarkonium itself in the lab- 
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Figure 3. Differential cross-section of 3/ijj pro- 
duction at CDF with theoretical predictions for 
color-singlet and color-octet model production [6] 



oratory frame. The polarization parameter a, de- 
fined as a = {<7t — 2ol)/((7t + 2ox), is equal to 
+1 for transversely polarized quarkonia produc- 
tion (helicity ±1). ot and <jl are the transverse 
and longitudinal cross sections. For longitudinal 
(helicity 0) polarization, a is equal to —1. Unpo- 
larized production consists of equal fractions of 
helicity states +1, and —1, and corresponds to 
a = 0. 

The previous Tevatron measurements were lim- 
ited to below around 20 GeV, where the polariza- 
tion is best predicted and the theory most under- 
stood. At ATLAS we aim to measure the polar- 
ization of directly produced prompt quarkonium 
in the pr region up to ~ 50 GeV and beyond 
with extended coverage in cos#*. This will allow 
for improved fidelity of efficiency, provide better 
discrimination of longitudinal and transverse po- 
larizations and therefore reduced systematic un- 
certainties. 

Evidently the polarization measurement can 
significantly suffer from low | cos 9* | acceptance, 
and hence from difficulties in separating detector 
efficiency corrections from polarization state ef- 



Figurc 4. Polarization of 3/tJj production as a 
function of pr at CDF, (data points), curves for 
limit cases of the fcr factorization model, and a 
band for NRQCD predictions (from [7]). 



fects. The di-muon trigger requirement for both 
muons to be above a certain pt threshold (/i6/i4) 
reduces to a minimum the acceptance at large val- 
ues of | cos 9* | , where the difference between var- 
ious polarization states is the more pronounced 
(see for example the two polarization states on 
the lower plots of Figure [5]). The acceptance of 
the single muon trigger sample (/xlO) is very dif- 
ferent. Here the efficiency is higher at large val- 
ues of | cos 0*| and drops around zero. Figure [5] 
(on the top) demonstrates how the two samples 
complement each other mainly at the low px re- 
gions, while at high pt the two triggers increas- 
ingly overlap, thus allowing for a cross-check of 
acceptance and efficiency corrections. The px dis- 
tributions of both samples, /x6/x4 and /^10 events, 
were appropriately combined. The combined dis- 
tribution <iiV raw jd cos 9* , was corrected according 
to equation [3l 

dN COT _ 1 dN raw 

dcos9* A(pr, cos 9*) ■ Ei ■ 62 dcos9* 

Here E\ stands for the trigger and reconstruc- 
tion efficiency, while e-i denotes the efficiency of 
background suppression cuts for each sample, and 
A(pt,cos9*) is the kinematic acceptance of the 
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trigger selection. The results obtained when fit- 
ting an unpolarized sample (a = 0) with statis- 
tics corresponding to 10 pb _1 are: a — 0.156 ± 
0.166,-0.006 ± 0.032,0.004 ± 0.029,-0.003 ± 
0.037, -0.039 ± 0.038 and 0.019 ± 0.057 corre- 
sponding to the pt bins as in Figure [5j The preci- 
sion of the cross-section derived from the normal- 
ization factor are ±4.35 in the first slice and de- 
creasing from ±0.09 to ±0.04 with increasing px 
slices. Repeating the same study with T we get 
back numbers that are consistent with (a = 0) 
but with larger errors running from ±0.17 to 
±0.22. To further check the ability to measure the 
polarization we reweighted our raw Monte Carlo 
(MC) distribution to emulate a transversely po- 
larized (a = +1) and a longitudinally polarized 
(a = —1) J/ip samples. The same analysis was re- 
peated and the results can be seen on the bottom 
of Figure [5j 

A complementary method for polarization 
measurement has been examined as well [5] . Here 
a linear combination of three MC generated sam- 
ples representing longitudinally polarized, trans- 
versely polarized and unpolarized events were fit- 
ted to the angular distribution of the data. The 
measured values of a in bins of px were well con- 
centrated around the initial value of the polariza- 
tion that was set in the MC-data sample. 

5. Analysis of \ production 

A sizeable fraction of prompt J/ip and T are ex- 
pected to originate from radiative decays of heav- 
ier states, \c and xt- These states have even C 
parity and therefore have a strong coupling to 
the color-singlet two gluon state. About 30 to 
40% of J/tp in our signal will come from decays 
of Xc — > J/ip + 7- Unfortunately, the energies of 
these photons tend to be quite small, and the abil- 
ity of ATLAS to detect these photons and resolve 
various x states is rather limited. For Xc recon- 
struction, each selected quarkonium candidate is 
combined with every reconstructed and identified 
photon candidate in the event, and the invari- 
ant mass of the /i/17 system is calculated. The 
H/i'j system is considered to be a x candidate, 
if the difference between the invariant masses of 
the fj,fj,"f and /i/i systems lies between 200 and 



800 MeV, and the cosine of the opening angle be- 
tween the J/ip and 7 is larger than 0.97. 

We fitted three guassians to the difference in 
invariant masses of the ///K7 and /in, measured in 
those x candidates. The MC input amplitudes of 
the peaks were reproduced reasonably well [5]. 
The overall Xc reconstruction efficiency is esti- 
mated to be about 4%. 

6. Physics reach with early ATLAS data 

The ATLAS B-Physics program will run from 
the earliest days and will pursue indirect searches 
for New Physics via B-hadrons decays. An ef- 
ficient, fast and clean di-muon trigger scheme 
will allow ATLAS to collect large samples of B- 
hadrons and quarkonia throughout the lifetime of 
the experiment. Reconstruction of exclusive B- 
meson decays with two muons in the final state, 
especially with J/ip, is foreseen already with the 
first pb^ 1 . Early B-Physics data will provide valu- 
able information on the detector performance and 
will allow calibration studies in support of New 
Physics searches. Meaningful quarkonium studies 
performed with early data are expected to have a 
sufficient reach to make authoritative statements 
about the underlying production mechanism and 
provide cross-section measurements in this new 
energy regime. Quarkonium spin- alignment mea- 
surements at ATLAS will have the capability to 
distinguish between various production models of 
quarkonium. Even in the early stages ATLAS is 
expected to provide competitive measurements of 
the J/ip and T polarization - px dependent val- 
ues. With several million J/ip fifj, decays, and 
better understanding of the detector, Xc — > Ji>l 
should become observable. 
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Figure 5. Six figures at the top show kinematic acceptances of the /x6/x4 (solid line) and /ilO (dashed 
lines) generator level cuts, calculated with respect to the sample with no generator level cuts on muon 
pr, in slices of J/t/j transverse momentum: left to right, top to bottom 9-12 12-13, 13-15, 15-17, 17-21, 
above 21 GeV. 

Six bottom figures show combined and corrected (equation [3]) distributions in polarization angle, for 



